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Abstract 

To investigate an effect of social interaction on the bystanders' intervention in emergency situ- 
ations a rescue model was introduced which includes the effects of the victim's acquaintance with 
bystanders and those among bystanders from a network perspective. This model reproduces the 
experimental result that the helping rate (success rate in our model) tends to decrease although 
the number of bystanders k increases. And the interaction among homogeneous bystanders results 
in the emergence of hubs in a helping network. For more realistic consideration it is assumed 
that the agents are located on a one-dimensional lattice (ring), then the randomness p S [0, 1] 
is introduced: the kp random bystanders are randomly chosen from a whole population and the 
k — kp near bystanders are chosen in the nearest order to the victim. We find that there appears 
another peak of the network density in the vicinity of A; = 9 and p = 0.3 due to the cooperative 
and competitive interaction between the near and random bystanders. 
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I. INTRODUCTION 



The concepts and methods of statistical physics and nonhnear dynamics are apphed 
to investigate the social, economic and psychological phenomena Among the 

interesting subjects that have attracted physicists are the o p inio n dy namics including 



voting process 



10, 



n 



12 



EI. 



Social impact theory 



P, y| and social impact theory 
stemmed from the bystander effect by which people are less likely to intervene in emergencies 
when others are present than when they are alone as a result of the inhibitory interaction 
among bystanders [3| • 

From the laboratory experiments about the emergency situations we can gain an insight 
into this effect. When tested alone, subjects behaved reasonably and the response rate was 
high. However the rate was significantly depressed when they were with other subjects. 
Subjects with others were unsure of what had happened or thought other people would or 
could do something. In another experiment subjects who were friends responded faster than 
those who were strangers. The subjects who had met the victim were significantly faster 
to report victim's distress than other subjects. And the degree of arousal that bystanders 
perceive is assumed to be a monotonic positive function of the perceived severity and clarity 
of the emergency, and bystanders' emotional involvement with the victim which is also 
to be considered in an abstract way in Section |n] 

In order to investigate the social phenomena as complex systems more precisely we adopt 
the network point of view by which it means that a social system consists of the interacting 
agents, where each node and link of the network represent an agent and a relation or inter- 
ac«o„ between a pa. of age.s respectively U^UU- A nu.be. of ptopetties about 
the real world networks such as social, technological and biological ones, have been revealed 
and investigated. Two of the main features of the real world networks are the small world 
effect and the high clustering to be considered in this paper by introducing a randomness 
p, the fraction of the randomly chosen bystanders to the k bystanders per accident, which 
plays a similar role to that of Watts-Strogatz model [2^ . The original dimensionless model 
for the bystander effect is extended to the more realistic and general one in Section lllll 
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II. RESCUE MODEL 



Recently in order to investigate an effect of social interaction on the bystanders' interven 



221. The model includes 



tion in emergency situations a rescue model (RM) was introduced 
the effects of the victim's acquaintance with bystanders and those among bystanders. The 
RM focuses on the relations between agents rather than on agents themselves, so defined is 
a relation spin between two agents as aij whose value is 1 if agent i (agent j) has succeeded 
in rescuing agent j (agent i) , and otherwise, is symmetric and can be interpreted as an 
element of adjacency matrix of helping network. Each agent i has its intervention threshold 
Ci over which that agent can try to intervene in an emergency situation. 

At each time step an accident happens which consists of the degree of the clarity or sever- 
ity of accident represented as a random number uniformly drawn from [0, 1], a randomly 
chosen victim v and k bystanders which are also randomly chosen from a population. Ny 
denotes the set of k bystanders. For each bystander i the degree of willingness to intervene 
Xvi is calculated: 

Xyi{t) = + aa^iit) + P ^ (2aij(t) - 1) - (1) 

Only one bystander i & with the largest value can intervene per accident, which can 
be called the intervener selection rule. If we assume that the response speed of bystander 
i is exponential in Xyi, the selection of the bystander with the largest Xm is justified. Ad- 
ditionally, once one bystander intervenes, the pressures on the others will disappear. Then 
the adjacency matrix is updated as following: 

a„(it + 1) = e(x„(t)) (2) 

where 9{x) is a heaviside step function. If Xm > 0, the rescue succeeds and then for the 
bystander i who intervened, the gets the new value of one. In case of x^ < the 
rescue fails and then the am gets the new value of zero, a represents the degree of victim's 
acquaintance with bystander, so can be called an acquaintance strength. The third term of 
Xyi is related to the interaction among bystanders. 2aij — 1 gives 1 if one bystander has 
succeeded in rescuing the other or —1 otherwise. There does not exist any neutral relation 
here. /? is used to tune the strength of coupling so can be called a coupling strength. Among 
them the main control parameter is the number of bystanders k. As observables we adopt 



the network density 17[ (helping rate in Ref. 2M) success rate respectively: 
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Sfc = ^5^^(x„(t)). (4) 

t=o 

In other words the success rate is defined as the number of successful interventions divided 



by the total number of interventions. Although the network density can be regarded as a 
kind of help ing rate, the success rate is closer to the helping rate defined in the experimental 
studies in a sense that the intervention may either succeed or fail without changing the 
network density. We fix Cj = c = 0.25 for all i according to the experimental result [ij] that 
70 ~ 75% of isolated subjects intervened and c does not change through this paper, which 
means we consider a population composed of homogeneous non-adaptive agents. Finally, 
the initial conditions are aij = for all pairs. 

At first let us consider the case without the coupling effect among bystanders, i.e. j3 = 0. 
Generally, an equation for the network density can be written as (22j | 

dak{t) 



dt 

where 



TVo^i - W^i^o, (5) 



H^o^i = (l-c)(l-afc(t))^ 
iyi_o = {c-a){l-{l-ak{t)f). 

Wq^i denotes the probability of creating a new link between the victim and the bystander 
and Wi^o does that of eliminating the existing link between them. The stationarity condition 
for ttfc yields 

which says is a monotonically decreasing function of k. In the numerical simulations ak{t) 
fiuctuates around since the links are added or removed with finite probabilities 1 — c and 
c — a respectively. As k increases, so does the probability that two connected agents, one 
as a victim and the other as a bystander, get involved in an accident again. According to 
the intervener selection rule one of the bystanders connected with the victim must intervene 
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and thus there is no reason for the increase in at according to k. Consequently the helping 
network gets sparse with the number of bystanders. 

An equivalent of the success rate defined in Eq. is given by 

= VTo^i - W^^o + 1 - (1 - auf = (7) 

1 — a 

where we used the stationary solution for in Eq. turns out to be independent of 

k and of the network density too. In fact, for the sparser network each link should bear the 
more burden on the intervention to ensure the success rate constant of k. From a viewpoint 
of the uncertainty of a victim's receiving help from the bystanders corresponds to the 
cost that the victim should pay to minimize the uncertainty. 

If the coupling effect among bystanders is taken into account, then from the definition 
of Xm the condition for the successful intervention can be obtained by a mean-field approx- 
imation, i.e. the substitution of for each aij-. 

Xvi = qv + OiCik + P{k - l)(2afc - 1) - c > 0, (8) 

or 

qv > -(« + mk - + I3{k-l) + c = ql. (9) 

At any time step, when given au the success rate corresponds to 1 — g*. In case with /? > 
there appear two transition points ki = + 1 and ^2 = + 1 (see Fig. Q). At k = ki 
for any accident the rescue succeeds, Sk = 1, if and only if ak = 1 while aX k = k2 for any 
accident the rescue fails, = 0, if and only if = 0. In the range of fci < A; < A;2, it 
is evident that ~ ql ~ for ajt > | + ^^2p{k-i) ■ O^'^^ ^ larger than q;/2, then 
the helping network is so dense that the probability that the bystander who has not been 
connected with the victim intervenes is extremely low, so is the possibility of creating a new 
link. One can expect that ak{t) increases since ~ 1, but very slowly since the network is 
sufficiently dense. 

Given Wi^q = we can calculate the time evolution of akit) by considering only the 
Wq^i. In case that the victim is not connected with any of bystanders, if we assume that 
at least one bystander is connected with all other bystanders, then for ki < k < k2, 

^ = Wo^, = {l-c-P{k- 1))(1 - ak{t)r 
at 

= (3{k2-k){l-ak{t)f. (10) 



Taking afc(t = 0) = as an initial condition yields 

ak{t) = 1 - [P{k2 ~k){k- l)t + . (11) 

Therefore ak{t oo) = 1 for fci < A; < k2- This solution represents the monotonically 
increasing behavior of the network density with time step and the k dependence as well. 

The time series of a^it) shown in Fig. |21 verify the above arguments except that the 
transition occurs at /c = 18 larger than ki expected by the mean-field approximation because 
of the finite size effect. One can see from the Fig. IHlthat as the system size increases, the 
transition point approaches ki = 16. Additionally afc(t) exhibits the punctuated equilibrium- 
type behaviors at k slightly smaller than the transition point ki, which will be revised in 
relation to the network viewpoint. 

Figure El shows the numerical results for and a^, both of which decrease until k reaches 
9 to 12. This tendency can be interpreted as the bystander effect in that the bystanders are 
less likely to intervene in emergencies (succeed in rescuing the victim in our model) when 
others are present than when they are alone. Contrary to the case with j3 = the decreasing 
ak according to k has an additional negative effect on the coupling among bystanders due 
to the positive (3, thus lowers the degrees of willingness Xyi in Eq. (jH)) and consequently Sk- 
However, Sk and are getting large as k approaches ki because of the excitatory coupling 
among bystanders. 

Next, let us focus on the effects of the acquaintance strength a and the coupling strength 
f3 on the structure of helping networks. If a = /3 = 0, since the degrees of willingness for all 
bystanders are the the helping network shows a completely random 

structure. If we consider the acquaintance effect, i.e. a > 0, the probability that two 
connected agents get involved in an accident again increases. Therefore a has a 'fixation' 
effect on the helping network. If the coupling among bystanders is taken into account, i.e. 
f3 > 0, the probability that the bystander connected with more other bystanders is more 
likely to intervene in an emergency, thus (3 has an effect of 'preferential attachment (PA)' on 
the helping network. As a result of the PA the heterogeneous hubs and hierarchical structures 
emerge from the homogeneous non-adaptive population. The PA has been investigated and 
summarized in Refs. [3, Q • 

In addition, interestingly the above punctuated equilibrium-type behaviors of afc(t) in Fig. 
121 accompany the rises and falls of hubs when they undergo the slow saturations punctuated 
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by the abrupt declines. The nontrivial total collapse of helping networks can result from 
the chain reaction between the effect of cutting links due to the rescue failure and that 
of the rescue failure due to the increasing negative interaction among bystanders. This 
phenomenon is very different from those of other cases in which once one agent becomes a 
hub, it lasts forever. 



III. RESCUE MODEL IN A SMALL WORLD 



In the previous section we ignored the spatial property of the system which does matter 
in realities. For more realistic consideration the randomness p G [0, 1] is introduced: when 
assumed that the agents are located on a one- dimensional periodic lattice (ring), the kr = kp 
among k bystanders are randomly chosen from a whole population, which can be called the 
random bystanders, and the kn = k — kr bystanders are chosen in the order in which they are 
nearest to the victim in the Euclidean space, which can be called the near bystanders. The 
near bystanders are to the local neighborhoods what the random ones are to the travelers 
from other places and so on. The randomness p makes the long-range interaction possible 
and plays the similar role in our model to the randomness defined as a control parameter of 
Watts-Strogatz small world networks 2l| . 



A. Agents in the One-dimensional World 

Let us first consider the case with p = which means that all the bystanders are the near 
ones. In case of even k, one half of bystanders are left to the victim and the other half are 
right to the victim. In case of odd k, k — 1 bystanders are chosen as for the case of even k 
except that the side of the last (farthest) bystander is chosen randomly, that is, left or right 
to the victim. We define a new observable as following: 

2/fcW = Tr/^E««^W (12) 

where [x] is a ceiling function and the denominator is the maximum number of links limited 
by the locality of interaction. By the definition of a^, 

2 2\k/2'\ 

^ ' i<j 



In case with P = 0, since equations ©-(0) for the case without locahty, i.e. for aA:,p=i, are 
vahd for i/k, it is natural to regard yt as ak,p=i. Thus for small values of k the network for 
the case limited by locality becomes much sparser than that for the case without locality. 
Interestingly Sk,p=o turns out to be independent of k again, precisely Sk,p=o = Sk,p=i = jz^- 
Similar to the reason for the k independence of Sk, in one-dimensional rescue model the 
probability that two connected agents get involved in an accident again is very high, thus the 
helping network gets sparse and as a result each link bears more burden on the intervention. 

In case with f3 > the helping networks in the one- dimensional world consist of a few 
hubs induced by the PA effect and their peripheries. The number of hubs amounts to about 
N/k and the number of peripheries per hub does to about k as shown in Fig. HJ Once the 
degrees of any agents become larger than those of others by chance, they eventually grow to 
the hubs and intervene in emergencies involved with their own peripheries and vice versa, 
which forms some kind of helping communities. In addition although the helping network in 
Fig. E] (b) does not show the scale- free behavior of degree distribution its backbone structure 

n n 

bears some resemblance to that of the structured scale-free network |22i] (see Fig. 8 in |2J] 
for comparison). 



B. Agents in the Small world 

The network densities and the success rates are scanned for the entire ranges of the 
number of bystanders k and the randomness p. When /3 = the numerical results depicted 
in Fig. El show the trivial behaviors. For each k, according to p the network density a^^p 
leaps from afc,p=o to about ak,p=i as soon as at least one random bystander appears, where 
kpc = 1 OT Pc = 1/k. For the values of p > Pc the network densities rarely change regardless 
of p, which implies that what is relevant is only whether the interaction is local or not and 
the other factors do not matter. The uncertainty of receiving help is maximized at k = 1 
and p = 1, where there is only one bystander chosen completely randomly per accident. 
Therefore the network density should be maximized to ensure the success rate. As seen in 
Fig. El(b), Sfc_p is independent of k as well as of p since the coupling effect among bystanders 
is not taken into account. 

If the coupling effect among bystanders is considered then an interesting phenomenon is 
observed in Fig. IHl that is, there appears another peak of ak^p and Sk,p in the vicinity of 
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k = 9 and p = 0.3. To understand the new peak in the intermediate range of both k and 
p we focus on a cooperative or competitive interaction between two groups; group of near 
bystanders and that of random ones. Based on the relevance of positive (3 we conjecture that 
the clustered structure of near bystanders is essential to enhance the random bystanders' 
intervention and their possibility of success so that the overall network density dominated 
by the nonlocal links can grow to a large value. When k fixed such as 9, for small p (large kn) 
local interactions among clustered near bystanders dominate nonlocal interactions among 
random ones and those across near and random ones, where the locality restrains the network 
density from getting large. For intermediate p and kn the random bystanders are benefited 
from the clustered near ones by making use of the existing links across near ones and random 
ones near to the near ones and then the range of interaction is expanded to the whole system 
after all. Consequently the network density becomes large. Finally, for large p (small kn) 
the near bystanders rarely cluster so that the random ones cannot be benefited from the 
clustering of near ones, hence the network density will converge to cifc,p=i rather than increase. 
To verify the above conjecture, at first two network densities are newly introduced: 



where d{i, j) gives the shorter distance on the ring between agent i and j and the superscripts 
n and r represent the near and random bystanders respectively. Nkn in Eqs. fll4|) - (fT3j) is the 
maximum number of links among near bystanders. From the extensive numerical simulations 
it is found that for the values of k and p other than the new peak region and the growth 
region of ki < k < k2 and p near to 1, both a^'^\t) and a'^'^\t) fluctuate around some 
values. On the other hand, for the new peak region, precisely at = 9 and p = 0.3, in 
Fig. [7| (a) a^"^(t) jumps to about 0.15 very quickly and fluctuates around that value for 
a while until a^^^t) grows exponentially to exceed a^"'\t), which is called the intersection 
point. a^'^\t) increases fast then saturates while a*^")(t) also increases a little. To figure 
out whether two bystander groups cooperate or compete we calculate the cross-correlations 
between two network densities before and after the intersection point. Before that point the 
cross-correlation is 0.5236 while after that point it is —0.3565. In the early stage the near 
bystanders help the random ones intervene but once the random ones dominate the near 




d{i,j)<ki 



(14) 




(15) 
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ones, two groups compete for the chance of intervention. 

For more specific investigation we calculate the averaged degree of willingness of random 
bystanders Xyr{t) as a function of time and that of near ones respectively as shown in 

Fig. [7|(b). Similar to the time evolution of a^'''\t), ^^.(t) increases from a lower value than 
Xyn{t) then exceeds it at the intersection point while x^nit) fluctuates around some value. 
Xyr larger than x^ implies the more chance of random bystanders' intervention than that 
of near ones' intervention. As a result the success rate of random ones (0.8531 ± 0.0101) is 
also slightly larger than that of near ones (0.8212 ± 0.0221). The difference of two averaged 
degrees of willingness affects which kind of bystanders are more likely to intervene and 
succeed, which affects the network density and the success rate successively. For general k 
by mean-field approximation the degrees of willingness are given by 



Xvr = qv + aak,p + (3{k - l){2ak,p - 1) - c, 
Xvn = qv + Oiak,p=i + (3{k - l){2ak,p=i - 1) - c, 

where we have assumed that the probability that a random bystander is connected with 
any other bystander is the same as a^^p and for the near bystanders the probability to be 
connected with any other one is Uk = a,k,p=i- From these the mean- field network density is 
obtained: 

ak,p -«^-P=^ + « + 2/?(fc-l)- ^^^^ 

The numerical results in Fig. [7| (c) for Eq. ()16p with k = 9 indicate that the mean-field 
approach works. 

For the other values of k what happens according to the randomness p? For smaller 
values of k and p > Pc the clustering of near bystanders rarely contributes to the random 
ones' intervention because the network is relatively dense to ensure the success rate for small 
k. Therefore it is the locality that is relevant for the results as for the case with /3 = 0. For 
larger values of k, especially larger than ki, we already know the network density goes to 1 
but very slowly when p = 1. Since for small p the network density converges to some value, 
there should be a transition line where the clustering cannot exist, i.e. kn < 1, equivalently 
p>l~l/k = p*. When p < p* and k > ki the locality limits the long-range excitatory 
interaction among bystanders and at the same time the probability that the victim and the 
bystander are already connected is relatively high due to large k so that it is difficult for near 
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bystanders to form an effective cluster and thus to enhance the random ones' intervention 
too. 

Next, as initial conditions we take the random networks whose network density a(0) varies 
from 0.05 to 0.5. Only for the small value of a(0), precisely 0.05, there appears the tiny peak 
in the intermediate range of k and p while no peak appears in the other cases with larger 
a(0). This is because the initial random structure does not allow any structural change for 
the clustering of near bystanders. In conclusion, the new peak in an intermediate range of 
k and p has been made sense by introducing the clustering effect of near bystanders. 

Finally it is also observed in Fig. IHlthat increasing a and /3 lead to the overall increase 
in the height of peak of network density then change the shape of peak from a sharp one 
to a plateau. Since a has a fixation effect on the existing victim-bystander link it cannot 
contribute to the creation of new link but only can lower the possibility of deleting the 
existing link. On the other hand (3 in relation to the bystander-bystander links can lead 
to the creation of new link between victim and bystander and the deletion of existing link 
depending on the network density. Therefore around the peak increasing (3 enhances the 
successful interventions of random bystanders based on the clustered near ones, then raises 
the network densities, which are also preserved by increasing a. 

IV. CONCLUSIONS 

In this paper we have studied not only the original rescue model, which was introduced 
in order to investigate an effect of social interaction on the bystanders' intervention in 
emergency situations, but also the rescue model on a small world. The bystander effect 
has been successfully reproduced from numerical simulations and explained by the mean- 
field approximation. In general both of the local interaction and the increasing k reduce 
the network density since the victim has more chance to get involved in the acquainted 
bystander. However, it is found that when the coupling effect among bystanders considered 
there appears another peak of a^ p and Sk,p in the vicinity of A; = 9 and p = 0.3 for some 
given parameters, which results from the enhancement of nonlocal interventions based on 
the clustering effect of near bystanders. 

The relation spins aij compose the helping networks. The coupling effect represented by 
positive (3 induces the emergence of hubs from a homogeneous non-adaptive population. In 
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the original rescue model the rises and falls of hubs have been observed near the transition 
point ki and in one-dimensional world the whole population is divided into a few helping 
communities, each of which consists of a hub and its peripheries. Although we could not 

find any real world helping networks to our knowledge, these results give us an insight into 
the dynamics of helping behavior and networks. 
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FIG. 1: The diagram for the explanation of the existence of two transition points ki and k2, where 
ki = + 1 and k2 = + 1 from the mean-field approximation Eq. @. Here a = 0.1, /? = 0.01 
and c = 0.25 are used. 
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FIG. 2: The numerical results of afc(t) for k = 9 (lower gray line), 16 (upper gray line), 17 (lower 
black line) and 18 (upper black line) respectively. ag{t) fluctuates around some value. aiQ{t) and 
an{t) repeat the slow saturations punctuated by the following abrupt declines. ais{t) increases 
monotonically but very slowly and finally approaches 1 as expected in Eq. t\ll^ . Here = 100, 
a = 0.1, P = 0.01 and c = 0.25 are used. 
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FIG. 3: The numerical results of the network densities and success rates for = 100 (circles), 
N = 500 (squares) and N = 1000 (plus signs) respectively. Each data point of Sk is obtained by 
averaging over last 5 x 10^ time steps of entire 1.5 x 10^ (5 x lO"^) time steps for N = 100 (for 
N = 500, 1000). And for ak the time span is the same as Sk but each point is averaged over every 
10^ time steps. For k > ki ak{t) increases monotonically as expected in Eq. (Illj) . however the 
points are obtained for finite time steps. Here a = 0.1, fi = 0.01 and c = 0.25 are used. 
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FIG. 4: The numerical results of the helping networks in the one-dimensional world for various k. 
The helping network for k = 11 is drawn in a circular style (a) and redrawn in (b). (c) and (d) 
are for the case with k = 31 and (e) and (f) are for the case with /c = 55 respectively. The hubs 
emerge from the homogeneous non-adaptive population. Here N = 100, a = 0.1, f3 = 0.01 and 
c = 0.25 are used and the networks have been produced with the Pajek software. 
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FIG. 5: (Color online). The numerical results of the network density and the success rate for 
1 < /c < 30 and < J? < 1 when /? = 0. Here = 100, a = 0.1 and c = 0.25 are used. 
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FIG. 6: (Color online). The numerical results of the network density and the success rate for 
1 < /c < 30 and < p < 1 when /? = 0.01. It is found that there appears another peak of a^^p and 
Sk^p in the vicinity of A; = 9 and p = 0.3. Here N = 100, a = 0.1 and c = 0.25 are used. 
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FIG. 7: The numerical results for verifying the importance of clustered near bystanders, (a) Time 
evolutions of two network densities introduced in Eqs. (|14() - H15|) : a(")(t) for near bystanders (black 
line) and a^'^\t) for random ones (gray line), (b) Time evolutions of the degrees of willingness; 
Xvn{t) for near bystanders (black line) and x^rit) for random ones (gray line). For (a) and (b) each 
point is averaged over 4000 time steps, (c) The verification of Eq. H16|) with A; = 9 by numerical 
simulations; x^n (crosses), x^r (plus signs), a^^^ (squares) and a^^p (circles). Each point is averaged 
over 50 realizations after saturated. Since there is no random bystanders for p < 1/k and there is 
no near ones for p = 1, in these ranges of p either or x^m cannot be defined. Here = 100, 
a = 0.1, (3 = 0.01 and c = 0.25 are used. 
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FIG. 8: (Color online). The numerical results of the landscape of network density for 1 < A; < 20 
(horizonal axis) and < p < 1 (vertical axis), changing from the sharp peak to a plateau according 
to a and (i: (a) a = 0.05, (i = 0.005, (b) a = 0.05, f3 = 0.02, (c) a = 0.2, [i = 0.005, and (d) 
a = 0.2, (3 = 0.02. Each point is averaged over 10 realizations. Here = 100 and c = 0.25 are 
used. 
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